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Abstract 

The conversion of ultracold atoms to molecules via a magnetic Feshbach resonance with a sinu¬ 
soidal modulation of the field is studied. Different practical realizations of this method in Bose 
atomic gases are analyzed. Our model incorporates many-body effects through an effective reduc¬ 
tion of the complete microscopic dynamics. Moreover, we simulate the experimental conditions 
corresponding to the preparation of the system as a thermal gas and as a condensate. Some of the 
experimental findings are clarified. The origin of the observed dependence of the production effi¬ 
ciency on the frequency, amplitude, and application time of the magnetic modulation is elucidated. 
Our results uncover also the role of the atomic density in the dynamics, specifically, in the observed 
saturation of the atom-molecule conversion process. 

PACS numbers: 03.75.Nt, 34.50.-s, 34.20.Cf 
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I. INTRODUCTION 


The optimization of the methods used for the production of ultracold molecules, in par¬ 
ticular, the minimization of losses, is an objective of current research 1|. Magnetic-field 
variations across a Feshbach resonance (FR) |^-4| and photoassociation 5 Tj are the tech¬ 
niques most frequently applied to molecule formation. In particular, linear magnetic ramps 


17] are standardly used in both Fermi and Bose gases. Additionally, Ramsey-like tech¬ 


niques, implemented via magnetic-field jumps close to a FR, have been applied to produce 


coherent superpositions of atoms and molecules 


a 


18 


19| . Here, we focus on an alternative 


method based on a sinusoidal modulation of the field in a magnetic FR. This scheme was first 
realized in experiments on thermal and condensed gases of 85 Rb 20]. In them, a resonant 


behavior of the efficiency was observed at a modulation frequency matching the molecular 
binding energy. The experiments revealed also the nontrivial dependence of the efficiency 
on the application time of the perturbation: the molecule population presented damped 
oscillations in the initial stage followed by a slow growth, which, eventually, saturated. A 
preliminary analysis, based on a two-body model with a two-channel configuration, showed 
that the observed features could not be understood in terms of a Rabi-like oscillation. A 
more complete description, in particular, the inclusion of different sources of decoherence, 
was required [20]. Theoretical and experimental work followed. The role of the system 


preparation was analyzed 


2l|. It was shown that, for a thermal gas, the initial distribution 


of atomic states leads to a quasi-continuous set of oscillatory components in the evolution 
of the molecular mode, and, consequently, to dephasing. The confinement in an optical 


lattice was also considered 


22], A variation of the modulation technique [23], applied to 


a mixture of different atomic species, uncovered molecular production at subharmonics of 
the resonance frequency. A similar effect was found for a single atomic species 24]. Despite 


the advances, there are still open questions; here, to go further in the clarification of the 
physical ground of the method, we extend former work on its general characteristics. We 
evaluate the differential role of the essential components of the dynamics in the emergence 
of the observed features. First, from an analytical description, we identify the characteris¬ 
tics rooted in the two-body Hamiltonian dynamics. Then, multi-atomic effects are assessed 
through a nonlinear configuration-interaction approach. In this framework, we analyze the 
effect of losses. Finally, we deal with the implications of the system preparation. Our ap- 
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proach is primarily applied to the experiments of Ref. [20]. The different time-scales and 
the eventual saturation of the transfer are explained as a combined effect of decoherence and 
nonlinearity of the (reduced) multi-atomic dynamics. Moreover, we uncover the origin of 
the resonances observed at subharmonics of the molecular binding energy in Refs. [23] and 
[24], A discussion of phase-space dependent restrictions to the efficiency, relevant to quite 
general contexts, is also presented. 

The outline of the paper is as follows. A first model, which incorporates basic two-body 
aspects of the dynamics, is presented in Sec. II. In Sec. Ill, an (effective) nonlinear Rabi 
model is applied to deal with the multi-atomic characteristics of the process. Molecular 
decay is included in this approach. The experimental preparation of the system is simulated 
in Sec. IV. Finally, some general conclusions are summarized in Sec. V. 


II. A TWO-BODY DESCRIPTION OF GENERAL DYNAMICAL CHARACTER¬ 
ISTICS 

We consider a Bose gas of ultracold atoms in a weakly confining harmonic trap. It is 
assumed that, as in the mentioned experiments, the interaction strength is controlled by 
a time-dependent magnetic held via a FR. In the practical arrangements, the held was 
modulated as 


B(t ) = B 0 + B m sm(ujt). (1) 

The average value B 0 was taken close to a FR position B r . Additionally, in order to work in 
a perturbative regime, the modulation amplitude B m was hxed at a sufficiently small value, 
B m <C \B 0 — B r \. The experiments of Ref. [20] were carried out on 85 Rb. In Ref. [23], a 
similar magnetic-modulation technique was applied to a mixture of 41 K and 87 Rb. Instead 
of a detailed characterization of each scenario, we will build a general framework where 
the elements responsible for the diverse observed features can be identified. A local-density 
approximation is applicable to the usual experimental conditions. Here, we concentrate on 
the simplest case of a uniform system: we will show that, in it, some of the experimental 
findings can be reproduced. Consistently, the atomic energy levels, which are slightly spaced 
because of the weak trapping, will be approximated as a quasi-continuum. 
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A. The undriven system 


Let us first consider the system without the magnetic modulation, (i.e., with B m = 0.) 
The starting point in standard approaches is a description of the coupled atom-molecule 
system in terms of bare states , i.e., of the eigenstates in the absence of coupling [1], The usual 
picture incorporates the (closed-channel) FR state \R), and coupled to it, the (entrance- 
channel) atom state | S). Correspondingly, the Hamiltonian reads Hi = H 0 + V, where H 0 
stands for the interaction-less system and V represents the coupling term. Thus, 

H 0 \S) = E s \S ), 0<E S <oo, 

H 0 |.R) = — 6b 0 l-R), e_B 0 > 0, 

(R\V\S) = v(E s ), (2) 

where Eg denotes the energy along the atomic set and e# 0 represents the binding energy 
of the state | R) at the mean magnetic held Bo. By now, we assume that the interaction 
strength v(Eg) hardly varies with Eg. Moreover, we consider that the set of eigenstates of 
the complete Hamiltonian Hi, i.e., the dressed states, consists of a bound eigenstate | M), 
(the molecular state), with energy —Em , and a continuum of states |A) with energy E A . 
Hence, 


Hi \A) = E A \A), 0 <E a < oo, 

Hi | M) = —Em | M ), Em > 0. (3) 

Both, | M) and |A), are given by superpositions of | R) and \S). For held values close to the 
resonance position, | M) can signihcantly diher from | R), and, in order to properly describe 
atom-molecule transfer processes, the dressed-state representation becomes necessary. This 
is the case of the considered situations: the resonant character of the modulation applied in 
the experiments refers specifically to the approximate matching of Em and fuv. 

B. The driven system 

In the practical setups [2(j|, the magnetic held was connected through a trapezoidal ramp, 
as that outlined in Fig. 1. The hrst linear segment slowly drives the held from B(t — tj), (far 
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from the FR), to B$. In the central plateau , the modulation described by Eq. (1) is applied. 
Finally, a slow linear ramp drives B{t) back to its original value. This arrangement demands 
the generalization of previous theoretical approaches. Accordingly, we develop a method 
which combines the following elements. First, an adiabatic approximation connecting the 
bare and dressed representations will be used in our simulation of the (slow) linear segments 
of cond, in the central region, 

on 



Figure 1: A diagram (with arbitrary scales) of the experimental magnetic-field ramp. 

The modulation is straightforwardly included in the bare-state picture: we assume that 
the held variation merely alters 6b 0 , according to e Bo + hCB 0 B m sin(ut), where hCe 0 = 
[ifli-Bo * s & iven by the difference between the magnetic moments of the involved states. 
(We will take h = 1.) Hence, the driven system is described by the Hamiltonian H 2 = 
H\ + C Bo B m sin (cut) \R) (i?|, which, in the basis of dressed states is written as 


H 2 = [~€ M + Vm sin(cnf)] | M) (M\ + 

[E a + rj A {E a ) sin(a;t)] \A) (A\ + 

[v eff (E A ) sin (ut) | M) (A\ + h.c.], (4) 

where rj M = C Bo Bm \{M\ R) | 2 , r) A (E A ) = C Bo B m \(A\ R) | 2 and 

v ef f(E A ) = C Bo B m (M I R) (R\ A) . (5) 
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The perturbation, apart from introducing a time-dependent coupling between the states 
| M) and |A), leads to a time variation in their eigenenergies. Furthermore, it can also 
couple different |A) states. However, that induced interaction between dressed atom states 
is neglected because of the small magnitude of the involved state-projections and the fact 
that the driving field and the atomic energy spacing are far from resonance. 

A comparison with the Landau-Zener (LZ) model 25], 26|, used to describe molecule 


production and dissociation induced by linear ramps, is of interest. In the LZ approach, 
the linear variation of the energy mismatch between the ( diabatic ) bare states leads to 
the eventual occurrence of a level crossing. In contrast, in our approach, the (sinusoidal) 
perturbation is not large enough for leading to crossings. In both models, the ramp starts 
and ends sufficiently far from the FR for the bare states to be approximate eigenstates at 
the extremes. As opposed to the standard solving method of the LZ model, our approach 
requires the switch to the dressed description. We will see that many-body effects can be 
included in our description in a form similar to that used in the counterpart generalization 


of the LZ model 


26 


221 


C. The secular dynamics 


We tackle the dynamics by applying the unitary transformation 
U(t) = e 4( iwt +^ cos Od)lW(^l+^ c ° s OdO)Ol] 


( 6 ) 


where l (> 0) is the nearest integer to the quotient between Em and c o. No restrictions on the 
perturbation frequency are assumed by now: the value of cu that optimizes the conversion 
will be determined in the following analysis. (Our study is not limited to the spectral region 
correspondin g to 1 = 1: we will also deal with frequencies quasi-resonant with fractional 


values of em 


23 


24].) The transformed Hamiltonian H 2 = H 2 U — iU'U reads 


H' 2 — {—Em + lw) | M) (M | + Ea |A) {A\ + 


v e ff sin (cut) 


\M)(A\ J2 H) k Jk(C/uj)e-^ k+l ^ t + h.c. 


k=— oo 


( 7 ) 


where Jfc(cc) denotes the ordinary Bessel functions [28] and ( = t/m ~ Va = 

C Bo B m (\(M\ R) | 2 — | (A| i?) | 2 ). This Hamiltonian provides a framework for simplifying the 
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description; specifically, an analytical coarse-grained picture of the dynamics can be derived 
through the averaging of H! 2 . The effective detuning — 6m + has been minimized by the 
choice of l as the nearest integer to e^/w. (Later on, the displacement of the resonances by 
thermal effects will be evaluated.) H' 2 has a series of interaction terms oscillating with fre¬ 
quencies (k + l±l)cu, k G Z. The secular dynamics is then governed by the components with 
k = —l- F1. The terms with k ^ — l =F 1, barely affect the system evolution: they present fast 
oscillations, which, because of the small magnitude of their amplitudes compared with the 
involved frequencies, can be averaged out to zero. Hence, in the optimum frequency range, 
namely, for lu ~ 6m, the dynamics is approximately described by the effective Hamiltonian 


| M) (M | + E a \A) (A| + 
'vf„(E A ) | M) </l| + h.c.l , 


where 5 ^ = lu — 6m, and 


( 8 ) 


Veff = [ J ~l~ l(C/w) + J-1+ l(C/w)] V eff (9) 

is a renormalized coupling constant which incorporates the magnetic modulation. This pic¬ 
ture provides the following clues to a preliminary understanding of some of the experimental 
results. 

i) Each value of / defines a resonance, —6m + lco = 0, around which, the dynamics is 

governed by a different effective Hamiltonian. The associated description corresponds to a 
discrete state coupled to a quasi-continuum, ft can be anticipated that, for each the 

transfer of population becomes more effective as the detuning decreases, which is consistent 
with the observed quasi-resonant character of the conversion process. 

ii) From the properties of the Bessel functions, it follows that the relevance of the inter¬ 

action term in the different H^ff decreases as / grows. (This is particularly evident through 
the approximation ~ valid for small \C/^\, as corresponds to the 

practical conditions.) These results agree with the findings of Ref. [23]. There, the pro¬ 
duction of molecules was observed to present different maxima at frequencies approximately 
given as u ~ 6m/ 1, (l > 1), the heights of the maxima decreasing with l. Similar features 
were found in Ref. [24], 
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iii) Some comments on the dependence of on the system parameters are in order. Let 
us illustrate the discussion by focusing on v^ff, which corresponds to the condition u ~ e^- 
In the perturbative regime, (Cs 0 i? m <C to), the effective coupling strength, given by Eqs. (5) 
and (9), can be approximated as 

e™ ~ j 0 (c/w)C Bo b,„ (m\ r) (R\ A) . (10) 

The factor Jq^/uj) does not incorporate significant variations with the system parameters: it 
does not appreciably differ from unity, as can be seen through the approximation Jq^/u) ~ 
1 — |(C/cu) 2 , valid for small arguments. Hence, v^ff is determined by B m and by Cb 0 , 
and, also, by the product of state projections (M\ R) (i?| A). In particular, depends 
on Bq, not only through Cb 0 , but also via (M\ R) (R\ A). These conclusions are consistent 
with the results presented in Ref. [29]. There, as the mean magnetic held Bq was varied, 
(and, in turn, Cb 0 was changed), the behavior of the coupling strength was observed to 
significantly depart from the mere linear dependence on CB 0 B. m . (For the terms with l > 1, 
the contribution of the factors that contain the Bessel functions can be nontrivial.) 

iv) No phenomenological approaches have been made in our identification of the mech¬ 
anism responsible for molecule production. The atom-molecule interaction induced by the 
magnetic driving has been traced in our first-principle scheme. However, the specific form 
of the effective Hamiltonians is strictly valid for small amplitudes and for frequencies in the 
quasi-resonant ranges. Outside those regimes, the applied approximation can break down, 
and, it is necessary to go back to H! 2 in Eq. (7) to have an appropriate description of the dy¬ 
namics. In fact, given the nontrivial dependence of H! 2 on u and B m , qualitatively different 
responses can be expected from significant variations in those parameters. 

The descriptions corresponding to the different effective Hamiltonians are parallel. Hence¬ 
forth, we will focus on the dominant term Bpjf, which will be simply denoted as H e ff. 
Moreover, we write = Em and v^ = h e //- Therefore, the two-body dynamics of our 
system will be assumed to be governed by 

H eff = E m \M) (M\ + E a \A) (A\ + 

[v e ff(E A ) \M) (A\ + h.c.]. (11) 

The quasicontinuum structure can be incorporated into the model by describing the 



atomic set in terms of the density of states and the coupling function. Note first that the 
possibility of varying uj allows a certain degree of control over the role of the continuum: 
different positions of | M) with respect to the atomic border, ( Em = u — Cm =3 0), can be 
realized by changing c o. For those different cases, the damping and energy shift of the discrete 
state are evaluated 30]. The results show that the features that we intend to explain, e.g., 


the different time scales or the saturation of the molecular population, do not seem to be 
rooted in damping due to the compact atomic-set structure. Therefore, we must explore 
alternative mechanisms: in the rest of the paper, we will evaluate the role of losses and 
many-body physics. Since we will neglect damping associated to the continuum structure, 
Eq. (11) will be regarded as defining an effective (linear) Rabi model involving | M) and 
an (isolated) generic state |A). This effective discretization implies the modification of the 
coupling constant to incorporate the characteristics of the density of states and the confining 
volume V. The specific form of the resulting dependence of the interaction strength on V 
is crucial to the role of the atomic density in the process. This issue was analyzed in Ref. 
[26] for a generic confinement. The results are applicable to our system: following them, 
we incorporate into Eq. (11) a modified coupling constant which depends on V -1 ^ 2 . [For 
simplicity, we maintain the notation v e ff(EA).\ In our reduced description of the atomic 
states, the variation of the energy Ea along the set is still contemplated. 


III. A NONLINEAR CONFIGURATION-INTERACTION APPROACH TO THE 
MULTI-ATOMIC DYNAMICS 

We turn to discuss how the above description can be generalized to study the multi-atomic 
scenario considered in the experiments. The generalization is necessary, in particular, for 
analyzing the dependence of the conversion efficiency on the atom density. 

A first-principle description of many-body effects is out of the scope of the study. Instead, 
we have opted for an effective reduction of the microscopic dynamics. From a qualitative 
picture, we can identify two elements that must be incorporated into the generalization of 
the two-body approach. First, we must address bosonic stimulation : since, for a condensate 
preparation, each of the N identical atoms can interact with N — 1 ~ N others to form the 
molecular species, a scaling of the interaction strength in the equations for the populations 
is pertinent. Second, the appropriate renormalization of the dynamical equations cannot 
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be simply implemented by introducing an N dependent amplification factor in the coupling 
constant. The correct procedure should account for the changing number of single atoms in 
the system as the conversion progresses. The inclusion of an effective interaction strength 
which depends on the atomic population implies the use of nonlinear equations, where 
unitarity in the evolution of the populations must be maintained. 

The methodology outlined above has been applied in previous studies. In early descrip¬ 
tions of transfers induced by linear ramps, the multi-atomic character was tackled by scaling 
the interaction strength in the equations obtained within a (standard) two-body LZ model. 


Specifically, in the expression 


C 


(as) 

M 


~ 1 — exp(— Wlz), which gives the asymptotic molec¬ 


ular population in the standard approach, the LZ parameter Wlz was replaced by Nwlz 


26]. [wlz grows with the inverse ramp speed. In the adiabatic limit, where Wlz —> oo, 


the model predicts a complete transfer of population. Furthermore, since Wlz scales also 
with the squared coupling matrix element, (and, in turn, with the inverse volume 1 /V), 
the effective parameter Nwlz grows linearly with the atomic density N/V. | In subsequent 
studies, it was also taken into account that the interaction strength changes as the molecular 


generation proceeds 


27]: a nonlinear configuration-interaction LZ model was implemented. 


Since the time scale for molecular decay is smaller than that of the complete conversion 
process, our model must also incorporate losses. Some numerical values illustrate this ar¬ 
gument: whereas the saturation of the process was found to occur on a scale larger than 
20 ms j20|, a typical time for molecular decay was estimated to be of the order of 1 ms 0- 
Loss mechanisms will be incorporated in a phenomenological way: we will describe them as 
a depletion of the molecular population characterized by the rate constant 7 f^ . 

Following the above lines, we assume that the state of the system, expressed as |\k(f)) = 
C A {t) 1-4) + Cm(£ ) | M), evolves according to a damped nonlinear Rabi model, built by 
extending the description given by Eq. (11), and explicitly defined by the set of equations 

H i 321: 


ijC A - Ul(N)C\C M 

( 12 ) 

i ^CM-\sr(N)c 2 A - 1 c M , 

(13) 


where A = Em — Ea = ca — (e« + Ea) and Q(N) = y/~Nv e ff. Since \v e ff\ 2 is proportional 
to the inverse volume, |fl(lV )| 2 scales with the atomic density. Further on, to deal with a 
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thermal gas, where no bosonic stimulation takes place, we will be back to the framework 
defined by Eq. (11). By now, in order to concentrate on the features specifically rooted in 
the preparation as a condensate, we consider the system in a pure dressed atomic state |A), 
i.e., we take C^( 0 ) = 1 , and calculate \CM{t)\ 2 - 


A. The undamped nonlinear Rabi oscillations 


For 7 = 0 , analytical solutions to the above set of equations can be obtained 


33]: the 


molecular population can be expressed in terms of the Jacobi elliptic function sn(a;, n) 281 
as 


\C M (t)\ 2 = pisn 2 (vWi |fI(lV)| t; rnj , (14) 

where p h2 = 1 + T J(l + i^yp) “ !, and m = PifP 2 - In Fig. 2, we compare 

the evolution obtained within our (nonlinear) scheme with that given by the standard Rabi 
model, which is recovered from Eqs. (12) and (13) by replacing the effective coupling 
constant f 1(N)C\ by Q(N). In contrast with the symmetric character of the (linear) Rabi 
oscillations, in the nonlinear case, the oscillations display asymmetry. Specifically, because 
of the incorporation of C^(i) into the effective interaction strength, the system evolves more 
slowly as the atomic population decreases. It spends more time near the molecular state, 
where \Ca\ = 0 . In particular, for A = 0 , and irrespective of the coupling-strength value, 
the evolution becomes aperiodic: once the system reaches t 
permanently. In this case, Eq. (14) is simplified to obtain 

\C M (t)\ 2 = tanh[|0(A)| f/2], (15) 

which reflects that, for Ea = Em, there is a complete and irreversible transfer of population, 
as opposed to the oscillations observed in the standard Rabi model. The time needed to 
complete the process is determined by the coupling strength: larger times are required as 
N\v e ff\ 2 decreases. This behavior is qualitatively similar to the ( irreversible ) evolution 
observed for linear ramps. Its persistence in the damped dynamics is the origin of some of 
the observed features. 

The dependence of the period and amplitude of the nonlinear oscillations on Q(N) and 


r_e molecular state, it stays there 
33] 


m 


ii 













A can be completely characterized from Eq. (14). Some general properties parallel those 
corresponding to the linear counterpart. Namely, as | A| grows, the amplitude and the period 
of the oscillations decrease. (This is illustrated in Figure 3.) Additionally, the amplitude 
in lensity. Depending on the 

m dynamics. 




Figure 2: Conversion efficiency \CM(f)\ 2 as a function of coupling time, as given by the nonlinear 
Rabi model [Eqs. (12) and (13) with 7 = 0] (solid line), and by the linear Rabi model (dashed 
line). The system is prepared in a pure |A) state and two values of the detuning are considered: 
A = 0 (a), and A = 1 kHz (b). (|D(A)| = 3 kHz). 
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Figure 3: Conversion efficiency |Cjvf(f )| 2 as a function of coupling time, as given by the nonlinear 
Rabi model [Eqs. (12) and (13) with 7 = 0]. The system is prepared in a pure |T) state and four 
values of the detuning are considered: A = 0.01 kHz (a), A = 0.5 kHz (b), A = 1 kHz (c), and 
A = 2 kHz (d). (|H(1V)| = 3 kHz). 


B. The role of molecular decay 


The occurrence of losses can significantly alter the conversion process. The total pop¬ 
ulation \CA{t )| 2 + |Cm (£)| 2 diminishes as the generated molecules decay and the emerging 
atoms escape from the trap. The magnitude monitored in the experiments was the dif¬ 
ference between the initial atomic population and that corresponding to a generic time, 
pj/(f) = 1 — \CAif) | 2 20 ], i.e., the sum of the molecular population and the fraction of 
atoms that have left the trap through molecular decay. In general, the form of P^{t) can 
be complex, with different stages and characteristic times. Still, some global properties can 
be identified. In a simple view, (which will be improved later on with the discussion of the 
phase-space proximity criterion), one can think that saturation occurs when all the atoms 
initially present in the trap have disappeared, and, consequently, when the maximum value 
P^iit) = 1 is attained. Since the atoms escape via transfer to the molecular state and 
subsequent decay, the scale for saturation depends crucially on the characteristic times of 
those component processes. Namely, when 7 is much larger than the effective conversion 
rate, the oscillations that could exist in the undamped system are absent, because, once 
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the molecules are formed, they suddenly decay. It is then the transfer that determines the 
scale for saturation. In the opposite limit, the evolution approaches the undamped oscil¬ 
lations: atom escape can be neglected at times much shorter than l/y. The behavior in 
a general regime incorporates traces of both extreme responses. This is apparent in Fig. 
4, where we depict P^{t) far from the above limits. The displayed features qualitatively 
reproduce the main experimental findings. A differentiated behavior is associated to A = 0, 
[Fig. 4(a)]: as no oscillations are present, |CA(t )| 2 continuously decreases irrespective of the 
relative magnitude of the characteristic times for losses and transfer, and, in turn, P^{t) 
monotonously grows. In contrast, as shown in Figs. 4(b) and 4(c), various stages appear for 
A^O. Damped oscillations, which are the remnants of the ones existent in the undamped 
model, can be seen if their period is smaller than I/ 7 . The frequency grows with both, 
the effective coupling strength and |A|. Because of the nonzero detuning, no inversion of 
population takes place; furthermore, for a sufficiently large |A|, the amplitude can be much 
smaller than 1 . Once the oscillations have decayed, the monotonous increase of P$(t) sets 
in, and the saturation is eventually reached. 

Since the appearance of the transitory oscillations demands an appropriate relative mag¬ 
nitude of the partial time scales, their stability can be affected by variations in the number 


of atoms in the different experimental runs 


29]. Actually, the changes in N affect VL(N), 


and, consequently, the Rabi period. The instability associated to thermal effects will be 
analyzed in the next section. 

Our picture applies specifically to the considered detection scheme. A modified approach 
can be necessary to deal with different arrangements. For instance, a continuous conversion 
process, which can be implemented via the replacement of the atoms that leave the trap, 
is described by adding a term of population gain to Eq. (12). The monitored magnitude 
P$(t) can then be identified with the molecular population, and the saturation corresponds 
to the balance between molecule formation and decay 0 


C. The phase-space proximity criterion 

In our model, as presented so far, P^{t) goes to unity, since, after a sufficiently long 
time, the initial atomic population has disappeared through molecule production and decay. 
However, in the experiments, the saturation was observed at a value smaller than 1. We can 
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think of having here a parallel of the reduced efficiency detected in linear-ramp transfers. 
As, in those processes, the efficiency was found to be determined basically by the density in 
the phase space, it was argued that a primary requirement for molecule formation could be 


the proximity of the atoms in position and momentum 


35]. This conjecture was validated 


by the reproduction of the experimental results through a simulation based on a stochastic 
phase-space sampling with an adjustable cut-off radius limiting the efficient region for as¬ 


sociation. (Additional support was given via coup 


36].) The wide applicability of this criterion 


37 


ed atom-molecule Boltzmann equations 


38] seems to convey the generality of 


the involved mechanisms. Still, the origin of the apparent phase-space depending charac¬ 
ter of the interaction is not clear 39]. From our approach, some clues to trace it can be 


given. The efficiency of the conversion depends ultimately on the magnitude of the matrix 
element v(E $) = (i2|V’|<S') (Eg = p 2 /m ; p denotes the relative linear momentum and m 
stands for the atom mass.) The variation of v(E$) with the atom-state characteristics is 
determined by the overlap between the wave functions ^r{t) = (r| R) and ip^(r) = (r| S), 
where r denotes the relative coordinate. (Since the operator V accounts for hyperfine in¬ 
teraction between the atom pairs in the open and in the closed channels, it affects only 
nuclear and electron spin variables, which are fixed in each channel.) The Franck-Condon 


can be evaluated under quite general assumptions 


40). It 


factor F(p) = f dr / il>f i (r)ijjg , \r) 
is found that F(p) becomes negligible for values of the momentum larger than pj j = h/d , 
where d represents the size of the bound state. Additionally, F(p) depends linearly on the 
inverse quantization volume, which implies, in both, condensed and thermal gases, a linear 
dependence on the density. For a condensate, the linear dependence on N/V emerges via the 
incorporation of the amplification factor N into the above expression of F(jp)\ for a thermal 
gas, the quantization volume in a two-body description is given by the inverse of the density 
ut- As the atoms are described in terms of free continuum states, (with phases modified 
by the scattering), the quantization volume determines the magnitude of the mean atomic 
distance (S\ r (S'). Hence, the requirement of having a significant Franck-Condon factor can 
be expressed as a phase-space proximity criterion: for the interaction to be relevant, the 
atoms must be in states with sufficiently small mean values of position and momentum. 
This supports the classical formulation of the criterion, namely, f.p < ~h, where, r and 
p, respectively denote the average relative distance and momentum, and x is a numerical 
factor obtained through the optimization of the simulation. From the quantum study, the 
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magnitude of x can be estimated via the expression (S\ r |/S') (S') ~ where |/S') rep¬ 

resents a state that effectively leads to molecular association. Taking (S'|p|S') ~ h/d and 
(S\r\S) d, the magnitude of the obtained value of \ corresponds to that of the classical 


estimate. 

The reduction in the efficiency imposed by the above criterion varies with the temperature 
and with the character, bosonic or fermionic, of the atomic sample. Here, given the variety of 
temperatures and condensate fractions realized in the experiments, we take into account the 
argument of phase-space limitations by considering that only a part of the atoms intervene 
in the conversion process |35|. As a result, our expression for the efficiency is rewritten as 
iVp(l — \CA{t)\ 2 )/N tl where N t denotes the initial number of atoms in the system, and N p 
corresponds to those that actually take part in the transfer process. In this phenomenological 
approach, applied to obtain the results of Fig. 4, we can identify the asymptotic efficiency 
N sion. (In the following, we 



Figure 4: Conversion efficiency pff (t) as a function of coupling time for the system prepared in 
pure | A) state for three different values of the detuning: A = 0 (a), A = 2 kHz (b), and A = 5 kHz 
(c). (|ff(lV)| = 3 kHz, 7=1 kHz.) [The inset shows the damped initial-stage oscillations for case 
(c).] (As in the following figures, the asymptotic efficiency has been arbitrarily fixed.) 
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IV. THE EFFECT OF THE THERMAL PREPARATION OF THE SYSTEM 


We evaluate now the dephasing effects that the system preparation can introduce in the 


dynamics. The majority of the experiments were done on uncondensed samples 


20|. In 


those arrangements, the starting point was a thermal distribution of atoms at a magnetic- 
field value far from the FR. At that distance from the resonance, the atom-molecule coupling 
is practically ineffective, and, consequently, the eigenstates of the complete system are ap¬ 
proximately given by the bare states for the initial held. The preparation is then a thermal 
mixture of \S) states. It can be assumed that, in the first segment of the trapezoidal ramp, 
(Fig. 1), each state \S) adiabatically follows the held, and converts, eventually, into its 
associated (B 0 ) dressed state, i.e., into one the states previously denoted as |A). Hence, the 
initial distribution is transformed into a mixture of |A) states. In the central region, the 
evolution of the atomic population is calculated applying the standard Rabi model. |The 
effective coupling constant in Eq. (12) is written as h e //(nr), since it depends on the atomic 
density tit through the quantization volume.] In the last segment of the ramp, the dressed 
states are adiabatically transformed into their related \S) states. 

From the propagation in the plateau of each component of the mixture, i.e., of each|A) 
state, we find |CA(t)| 2 . Then, the total atomic population is obtained by averaging |CA(t)| 2 
over the Boltzmann distribution of the original mixture of | S) states. The average can be 
simplified through the following lines. From previous characterizations of the coupling of a 
discrete level to a quasi-continuum, one can assume that the energy spacing between consec¬ 
utive states | A) is the same as that existent between their bare counterparts | S). Hence, an 
appropriate shift in the origin leads the energy E& of a (generic) |A) state to match the rela¬ 
tive kinetic energy p 2 /m of the pair of atoms in the original state | S). (Free-space conditions 
are being used to evaluate the energies and the density of states; p corresponds then to a 
continuous variable.) Namely, we can write E&(p) = p 2 /m, and, in turn, A (p) = Em~P 2 /m. 
Hence, the initial weight function simply corresponds to the Boltzmann distribution for the 
|A) states with shifted energies. (To indicate that the atomic population |CA(t)| 2 depends 
on p through A (p), we rewrite it as |CA(p, t)| 2 .) Since the statistical average is not affected 
by the displacement of the energy origin, the total atomic population can be approximated 


as 
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pJaW = /|c^ t )l 2 «P(-V/^ ( 16 ) 

where (3 = (/c^T) -1 . (In the parallel equation for the molecular population, the linear depen¬ 
dence on the atomic density, found in Ref [21] in the short-time limit, becomes apparent.) 
One can see that the averaging introduces decoherence: the (partial) oscillations present 
in \Ca(p, t) | 2 , rooted in nonzero values of A (p), dephase because of the almost continuous 
variation of Ea{p ) along the distribution. Then, in addition to damping generated by losses, 
there is decoherence due to the preparation of the system, its characteristic time being deter¬ 
mined by the width of the distribution, and, therefore, by the temperature. It is also evident 
that the initial-stage oscillations decay faster as T increases. Moreover, as T is changed, 
the distribution of detunings is altered, and, as a consequence, the average frequency can 
be shifted. The time required to reach the saturation regime depends, not only on intrinsic 
characteristics of the system, but also, on the width of the distribution. These results are il¬ 
lustrated in Figs. 5 and 6, where we represent the efficiency P^{t) = 1 — p^\(t) for different 
values of the threshold detuning A (p = 0) = Em = cu — 6m■ (From the experimental data, 
no information can be extracted on the volume and densities, and, therefore, on the effective 
coupling constants in the condensed and thermal gases; we have used arbitrary values for 
th 1 1 ' ’ “ ' loreover, equal phase-space 

re msidered temperatures.) 



Coupling Time (ms) 
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Figure 5: Conversion efficiency P^(t) as a function of coupling time for the system prepared as 
a thermal statistical mixture (T = 20 mk) for three values of the threshold detuning: A (p = 0) = 
5 kHz (a), A (p = 0) = 0 (b), and A (p = 0) = —2 kHz (c).) ( v e ff{n,T ) = 1-5 kHz, 7 = 1 kHz.) [The 



Figure 6 : Same caption as Fig. 5 with T = 40 mk 

Results for the association spectrum, i.e., for the efficiency as a function of the modu¬ 
lation frequency are presented in Ref. [20]. Some of the observed characteristics, like the 
changing form of the curves at different coupling times or the symmetry properties, must be 
explained. The multi-factorial character of the dependence of the efficiency on lu is central 
to the analysis, c 0 enters the populations through the coupling strength v e ff- The detuning, 
and, therefore, c 0 , affects also the amplitude of the oscillatory components of the evolu¬ 
tion. Moreover, the weights of the statistical mixture effectively modify the amplitudes. 
To isolate the thermal effects, we have worked with a frequency independent v e ff- (In the 
considered perturbative regime, this is a good approximation, since the factor J 0 ((/co) does 
not appreciably differ from unity.) In Figs. 7, 8, and 9, we represent the efficiency as a 
function of A = Em — Ea for a pure state, and, of A {p = 0) = Em for a thermal gas. The 
magnitude of the used parameters corresponds to the experimental conditions. Apart from 
common general features, the spectra present differential aspects specific to the preparation 
characteristics: 
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i) The resonant behavior observed in the experiments is apparent. For a condensate, 
(Fig. 7), the resonance occurs at the binding energy. For a thermal gas, (Figs. 8 and 
9), as the temperature grows, the resonance is shifted to increasing values of A (p = 0). 
Actually, in the condensate, the evolution corresponds to just one oscillatory term, which 
reaches its maximum amplitude when A = 0. In contrast, for a thermal mixture, there is a 
quasi-continuous set of oscillatory components in the molecular population. The energy of 
each state, and, in turn, the corresponding detuning, A(p), vary continuously along the set. 
Additionally, the statistical weights of the different terms change with the temperature: as 
T increases, the region of dominant contribution to the integrand in Eq. (16) is shifted to 
higher E&. As a consequence, the maximum of the transfer efficiency, which takes place for 
Em approximately matching the dominant region of Ea, occurs for higher values of E M , or, 
equivalently, for growing detunings. 

ii) For the preparation as a condensate, the spectrum is symmetric with respect to the 
position of the maximum (A = 0), since it is the magnitude of the (single) detuning that 
affects the efficiency. In contrast, there is no symmetry in the thermal case: because of the 
form of the distribution, the weights of the oscillating terms equidistant from the energy 
of the maximum are different. One can also understand that the decrease is faster in the 
region of negative threshold detunings, i.e., for effective molecular energies below the atomic 
set. In those cases, a zero value of A(p) = Em — p 2 /m cannot be reached. Then, there is 
no component in the mixture attaining its maximum (zero-detuning) amplitude. Moreover, 
since the magnitude of A(p) monotonously grows along the set, the amplitude of the corre¬ 
sponding oscillatory terms steadily diminishes. Consequently, the output of the transfer is 
less efficient in that region than in that of positive values of A(p = 0). 

iii) The observed time broadening of the spectra is also reproduced: as found in the 
experiments, the peak width increases with time and eventually saturates. Indeed, as the 
magnitude of the detuning increases, the conversion becomes less efficient, longer times being 
needed for molecule formation to take place. Once the asymptotic regime is reached, with 
all the available atoms converted into molecules, and, afterwards, decaying, the efficiency 
does not longer increase with time. 

Our results reproduce the form of the_experimental spectra 


predictions of previous theoretical work 


20 ] and coincide with the 


2lJ. 
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Figure 7: Conversion efficiency (t) as a function of detuning A/(27r) for the system prepared 
in pure |H) state for three coupling times: tf = 4 ms (a), tf = 18 ms (b), and tf = 34.4 ms (c). 



Detuning A(p=0)/(27t) (kHz) 

Figure 8 : Conversion efficiency Pff(t) as a function of threshold detuning A (p = 0)/(27t) for the 
system prepared as a thermal statistical mixture (T = 20 mk) for three coupling times: tf = 4 ms 
(a), tf = 18. ms (b), and tf = 34.4 ms (c). ( v e ff(nT ) = 0.9 kHz, 7 = 1 kHz.) 
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Detuning A(p=0)/(27t) (kHz) 

Figure 9: Same caption as Fig. 8 with T = 40 ink 


V. CONCLUDING REMARKS 

Our analysis opens the possibility of controlling the response of the system to the applied 
held. The partial analytical character of the results facilitates the design of strategies for 
improving the molecule production. The effective interaction strength can be regarded as a 
parameter of control. Appropriate modifications of it, through variations of the modulation 
characteristics, can be implemented to optimize the efficiency or to reduce the transfer time. 
One can also think of using a different confinement, (e.g., trapping in an optical lattice 
instead of the original harmonic potential), to have an additional element of control to vary 
the coupling matrix element and to reduce or suppress inelastic collisions. 

The analysis, although focused on a particular method of ultracold molecule production, 
can have general implications as it deals with fundamental mechanisms. In this sense, the 
identification of the role of the different components of the dynamics can give some clues 
to the description of related systems. The intended operative character of the study has 
implied simplifications in the presentation and tracing of our approach. We project a more 
complete account of the derivation of the model from a microscopic description. Another 
objective of future work is the evaluation in different contexts of our conjectures on the basic 


22 




mechanisms that can limit the efficiency. 
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